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Magnetization behavior of long nanowires witii square cross-section under influence of strong 
perpendicular spin-polarized current is studied theoretically. The study is based on Landau-Lifshitz- 
Slonczewski phenomenology. In the no current case the wire is magnetized uniformly along its 
axis. For small currents the wire magnetization remains uniform but it inclines with respect to the 
wire axis within the plane perpendicular to the current direction. With the current increasing the 
inclination angle increases up to the maximum value 7r/4. Further current increase leads either to 
saturation or to stable periodic multidomain structure depending on the wire thickness. For thick 
wires a hysteresis is observed in the saturation process under the action of current. All critical 
parameters of the current induced magnetization behavior are determined theoretically. The study 
is carried out both analytically and using micromagnetic simulations. 
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I. INTRODUCTION 

The usage of spin-polarized current is a convenient way 
to control the magnetic structures in nanowires with- 
out applying the external magnetic fields, that enables 
increasing of density of arrays of nanoscale elements in 
purely current controlled devices^'^. Usually, the current 
is passed along the magnetic wire, this is so-called CIP 
(current-in-plane) configuration. In this case the influ- 
ence of the current on the dynamics of domain walls is 
widely studied both theoretically and experimentally, see 
reviews [3-6]. 

In the last few years an interest to CPP (current per- 
pendicular to the plane) configuration of nanowires ap- 
pears, e.g. it was shown theoretically^ that in CPP case 
the velocity of the domain wall can be much higher than 
those for a CIP stripe with similar applied current den- 
sities. The similar domain wall motion was recently ob- 
served experimentally^'^ characteristic currents densities 
are much smaller than the one commonly used in CIP 
configurations. It was also shown^^di ^i^^t the domain 
wall can move even faster for the time alternating current 
in CPP wire structure. 

In this work we also consider the CPP nanowire con- 
figuration, see Fig. 1. We study possible stationary 
states and their stability under the current action. We 
found that there exists a critical transverse size of 
the nanowire. For sizes h < he & stable periodic do- 
main structure appears in the pre-saturated regime, for 
thicker wires h > he the process of transverse saturation 
by current has a hysteresis, i.e. the critical current of 
transition to the saturated state (in the process of the 
current increasing) is larger than the critical current of 
the saturated state breaking (in the process of the cur- 
rent decreasing). Recently we reported on the formation 
of stable periodic square vortex-antivortex lattice (vortex 
crystal) in the pre-saturated regime of the magnetic films 



of CPP configuration^^'^'^. The idea of this paper is to 
demonstrate that strong restriction of film size in one di- 
mension significantly affects the current-induced behav- 
ior of magnetization. In addition to the above-mentioned 
results we found that only stationary states take place in 
the nanowire with CPP, contrary to films where the es- 
sentially dynamic regimes, such as fluid-like and gas-like 
dynamics of vortices and antivortices appear for current 
decreasing^^'^'^. 

The paper is organized in the following way: in Sec- 
tion II we introduce a mathematical model of the one- 
dimensional magnetic nanowire which is based on the 
Landau-Lifshitz-Slonczewski equation. Here we adapt 
the spin operators representations of Holstein-Primakoff 
and Tyablikov for our classical system. Considering a 
case of soft ferromagnet we take into account exchange 
and dipole-dipole interactions and derive the correspond- 
ing Hamiltonian in the wave-vector space. In Section III 
we obtain two possible uniform stationary solutions and 
analyse their stability. The nonuniform stationary solu- 
tion - periodic domain structure, which arises as a re- 
sult of instability of the saturates state, we study in the 
Section IV. In whole, the magnetization behavior of the 
nanowire under the current influence is summarized in 
the Conclusions and mathematical details of deriving of 
an equation of motion and Hamiltonian are placed in two 
Appendixes. 

II. ONE-DIMENSIONAL DISCRETE MODEL 

We consider here a narrow nanowire of square cross- 
section whose transverse size h is small enough to ensure 
one-dimensionality of the magnetization, see Fig. 1. In 
other words we assume that the magnetization is vary- 
ing only along the wire - along x-axis. This assumption 
works for the cases when h is comparable or smaller than 
the characteristic magnetic length of the material, see 
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which can be index (coordinate) dependent. For this pur- 
pose, we also use the eigen-representation proposed by 
Tyabhkov^^ and finally, combining these two approaches, 
we obtain the following representation for the magneti- 
zation 



n^-4P^ + A„V^„V2-|^nP + c.c., (4a) 



FIG. 1. Geometry and notations of the problem. A long 
nanowire of square cross-section and thickness h. The spin- 
polarized current J flows perpendicularly to the wire. 



bellow Section II B. Total length of the wire L ^ h. The 
frame of reference is chosen as shown in the Fig. 1. The 
magnetic media is modeled as a discrete cubic lattice of 
magnetic moments Mi,, where v = a{vx, Vy, Vz) is three 
dimensional index with a being the lattice constant and 
^"2;, t'y, S Z. It is convenient to introduce the following 
notations: J^x ~ L/a is the total number of lattice nodes 
along X-axis and = 1? ja^ is the number of nodes 
within the cross-section square. 

We base our study on the one-dimensional discrete 
Landau-Lifshitz-Slonczewski equation: ^^^^^ 



where vector An is orthogonal to 7^ and has the following 
form 



=^ (cos e„ cos $„ -I- i sin $„, 

cos 0„ sin $„ — i cos $„ , — sin 9„) 



(4b) 



Using the representation (4) we obtain the following 
expressions for the magnetization components 



ml + iml = e^*"<{ sine„ (l - |V'„P) 
8 



(5a) 



/ ^ " /I I * ■ 2 ®n 

ipn COS -y-'^l^n sm — 



m„ = m„ X d£./dmn - je„m„ x [m„ x z], (1) ^ cose„(l - |i/;„n 



sin 0, 



which describes the magnetization dynamics under in- 
fluence of spin-polarized current which flows perpendicu- 
larly to the wire along f-axis as it is shown in the Fig. 1. 
It is also assumed that the current flow and its spin- 
polarization are of the same direction in (1). The index 
n ~ avx numerates the normalized magnetic moments 
m„ = Mn/\Mn\ along the wire axis. The overdot in- 
dicates derivative with respect to the rescaled time in 
units of (47r7oMs)~^, 70 is gyromagnetic ratio, Mg is the 
saturation magnetization, and £ = E/{4:7rM'^a'^l^s) is di- 
mensionless magnetic energy. The normalized electrical 
current density j = J/Jq, where Jg = 4TrM^\e\h/H with e 
being electron charge and h being Planck constant. The 
spin-transfer torque efficiency function e„ has the form 



(5b) 



It should be noted that for the case of absence of devi- 
ations we have ipn = and expressions (5) result in the 
magnetization orientated along vector 7„. The similar 
procedure we used in Ref. 13 for the case of the uniform 
stationary state which is oriented along the current. 

Substituting (5) into (1) enables one to proceed from 
the set of equations for the magnetization components to 
a single equation for complex valued function 



(6) 



(l + A-2) + (l-A-2)r 



(2) 



where the spin-torque term 3^^ l^as the following form 



where 77 is the degree of spin polarization and parame- 
ter A ^ 1 describes the mismatch between spacer and 
ferromagnet resistance^^'^^. 



A. Holstein-Primakoff-Tyablikov representation 

For the future analysis it is convenient to introduce 
complex amplitude ipn of deviation from the station- 
ary distribution of the magnetization (possible stationary 
states are considered bellow). For this purpose, we use 
classical analogue of Holstein-Primakoff representation^® 
for spin operators. Here we consider general case of an 
arbitrary form of the stationary state, where the magne- 
tization distribution is determined by the unit vector 



3^n — 



COSe„?A„ ( 1 - - \lpn\'^ 



- sin O, 



(7) 



7„ = (sin9„cos$„, sin 8„ sin <!>„, cos 8^ 



(3) 



Here the efficiency function e„ depends on accord- 
ing to (2) and (5b). For details of deriving (6) see Ap- 
pendix A. 

Since the Eq. (6) describes the deviation from a sta- 
tionary solution, it has a very convenient form for anal- 
ysis of stability of the given stationary state. Although 
the Eq. (6) can be used for arbitrary stationary state, in 
what follows we restrict ourselves to the case of spatially 
uniform stationary states, so that 8„ = 8 = const and 
$„ = $ = const. Here we study the linear stability, thus 
it is enough to use the linearized form of Eq. (6). It is also 
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convenient to proceed to the wave-vector representation 
using the discrete Fourier transform 



^•^^^ k 



with the orthogonality condition 



(8a) 
(8b) 

(8c) 



where k — is two-dimensional discrete wave vector, 
I £ Z, and A(fc) is the Kronecker delta. Applying (8) to 
the linearized equation (6) one gets the following equa- 
tion of motion in the wave-vector space 



9£0 



i]e 



sin 9 + cos Qij^k 



+ -(l-A-2) sin2e(V3fc+V>J 



(9) 



where e" = r// [(1 + A^^) + (1 - A^^) cos 6] is the spin- 
transfer torque efficiency function (2) for the case if^n = 
and EP denotes the harmonic part of the normalized 
energy. 



B. Hamiltonian of the system 



representation (5) to (10) and applying the Fourier trans- 
form (8) one obtains that the harmonic approximation of 
the normalized exchange energy reads 



(13) 



where the exchange length I — y^§'^3/i'2nM^a) deter- 
mines the scale of magnetization inhomogeneities. Here 
we neglected possible surface effects which can arise due 
to the different number of the nearest neighbors at the 
surface. The derivation of (13) is analogous to one pre- 
sented in Appendix Al of the Ref. 13. 

Let us proceed now to the dipole-dipole contribution. 
Using the fact that the magnetization depends only on 
x-coordinate one can present the expression (11) in form 



Mfd 



2„6 



E 



E 

S—x,y,z 



E 



3^, 



(14) 



Here we introduced the notations P = {vr^, Vy) and 
A = (m^i My) for the sake of simplicity. It is easy to 
show^" that -Bi/^^i^ = for the case of rectangular wire 
cross-section, the same is true for the cross-sections in 
form of disk (cylindrical wire) or ring (tube). Substitut- 
ing (5) into (14) and applying the Fourier transform (8) 
we obtain the following harmonic approximation for the 
normalized dipole-dipole energy 



We consider here the case of soft ferromagnet, therefore 
we take into account only two contributions into the total 
energy: E ^ E^x + Ed- Here 



Eex = -§^3 E ' "^^+6 



(10) 



is the exchange contribution, where 8 = a{5x, Sy, dz) is 
a three dimensional index which numerates the nearest 
neighbors of an atom, value of spin is denoted with § and 
3 > is exchange integral between two nearest atoms. 
The other term is the dipole-dipole energy 



Ea 



Mid 



2„6 



E 



{m^-m^) _ g (m 



(11) 



where we introduce the notation 

Hereinafter we will be interested only in the harmonic 
approximation EP of the normalized energy which in- 
cludes terms not higher than 0(17/)^). Thus 



(12) 



where and fi^ ^'"^ harmonic parts of exchange and 
dipole-dipole energies respectively. Substituting now the 



£2 = -^13^23^ sin eg(0) cos $(cos 9 cos $ + i sin $)^fc 



+ 29{^){'^'^^ Q cos <i> -I- i sin ^)'^ipk'4'-k 

+ I- [g{k) + 2g{0)] (1 - 3sin2 6 cos^ $)|7/'fcp + c.c. 



(15) 



where the function g{k) can be presented approximately 
as 



kh\ f kh 



(16) 



with Ii [x) and Ki [x] being modified Bessel functions of 
the first and second types perspectively. For the exact 
form of g{k) and other details see Appendix B. 



III. STATIONARY STATES OF THE SYSTEM 

The stationary magnetization distribution is deter- 
mined by Eq. (9) with ip = 



sm B = — = 0. 

V2 



(17) 
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FIG. 2. Components of total magnetization of the nanowire 
as functions of the apphed current. Panels a) and b) cor- 
respond to square permalloy nanowires with transverse sizes 
h = 5 nm and h = 15 nm respectively. The results of mi- 
cromagnetic simulations are shown by solid lines: thick and 
thin lines corresponds to simulations with periodic bound- 
ary conditions along x-axis and to the wire of finite length 
L = 1 jim respectively. The analytical solution (21c) is shown 
by the dashed line. The current density x is normalized in 
accordance with (20b). In region 1 the wire is magnetized uni- 
formly within the plane perpendicular to the current flow and 
inclined at an angle $ to the wire axis, see (20) and Fig. 3b. 
The region II corresponds to a periodic domain structure, see 
Fig. 4, and region III corresponds to a saturated state, see 
Fig. 3c. For the thick wires (panel b) the hysteresis region 
I-III arises, where the magnetization state coincides with the 
former state: either with the state I for the process of current 
increasing, or with the state III for the current decreasing 
from the saturation value. The hysteresis loop is shown only 
for component in order to not overload the figure. 



Substitution (15) into (17) results in the equation for 
ground states of the system for a given value of the cur- 
rent j 



sin 9 



i]e 



cos $ 



(cos Q cos $ — i sin $) 



= 0, (18) 



where we took into account that 17(0) = 1/12. 

Equation (18) has two solutions. The first one, namely 



sin 9 = 0, 



(19) 



describes a saturated state when the wire is uniformly 
magnetized along £-axis, i.e. along the current direction. 



The second solution of Eq. (18) reads 

cos 9 = 0, 
sin $ cos $ = — >f, 
2r; 



1 + A- 



(20a) 
(20b) 



In terms of magnetization components (5) the solution 
(20) has the following form 



0, 



■m^{x) = sin $(>f) 



1 + Vl - 4>f2 



^/l - 4x2 



(21a) 
(21b) 

-X. (21c) 



The solution (21) as well as (20) exists for the current 
interval < x < 1/2, that corresponds to varying of the 
angle $ in interval — 7r/4 < $ < 0. 

In the following we consider stability of each of the 
solutions (19) and (20). Let us start from the stability 
analysis of the solution (20). To obtain the equation of 
motion linearized in vicinity of the stationary solution we 
substitute (20) into (9) and (15), that results in 

- ^Vi. = ^ + ^ 4 r^A=2 (^^ + r-u) , (22a) 
£" = E |9(fc) sin^ ^x) (^fc^.fc + i^ir-k) (22b) 

k 

+ \M^{ek' + [g{k) + 2g(0)] [2 - 3 sin^ $(x)] } , 

where function sin^ <^{k) is determined by Eq. (21c). The 
linear Eq. (22a) has the solution 7/1^ = \I'-|-e^±^'^^*, where 
the rate function z± (k) is 



z±ik) 



1 r 
2 



r2 — 4a+a 



where 



a±(/fc) =£2^2 ^ [^(^j.) _^ 25(0)] [2 - 3sin2 $(; 

±3g{k)sm^^{>t), 
_ _>f 1- A-2 
^~2 1 + A-2- 



(23a) 



(23b) 



(23c) 



In accordance with Eq. (23a) the normalized current 5i 
plays a role of an effective damping2i, that is why the 
natural damping can be omitted in the original equation 
(1). When 



a+a < 



(24) 



the rate z± becomes positive what results in instability 
of the stationary solution (20). Analysis of (23b) shows 
that for h/£ < Gq with Cg ~ 17.37 the instability con- 
dition (24) is equivalent to the condition sin^ $ > 1/2. 
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Thus for the considered interval — 7r/4 < $ < or, in 
other words, for the current interval < >f < 1/2 the 
stationary solution (20) is stable. We do not analyze the 
case of very thick nanowires because the assumption of 
one-dimensionality of the magnetization does not work 
for the wire thickness hji^ 1. 

Thus, when the current is absent (>f = 0) the wire is 
magnetized uniformly along its axis. Adiabatically slow 
increasing of the spin-current from the value >f = to 
value >f = 1/2 leads to the homogenous inclination of 
the wire magnetization by the angle $(x) which accord- 
ingly changes continuously from value = to value 
<I> = — 7r/4, where the function $(>f) is determined by 
(20). The rotation takes place within the plane perpen- 
dicular to the current direction. The described regime 
takes place in region "I" and also in the hysteresis region 
"I-III" if the current is increasing, see Fig. 2 and Fig. 3. 
As one can see from the mentioned figures the analyt- 
ically obtained behavior of the magnetization is in full 
agreement with micromagnetic simulations^^. 

The critical current Kc ~ 1/2 in physical units reads 

Jc = ' (1 + A'') ■ (25) 

Thus dependence of the critical current Jc on the wire 
thickness h is linear one with the slope dependent on the 
parameter A. The corresponding dependencies Jc{h) are 
shown in the Fig. 3a by lines 1, 1', and 1" for different 
A-s. 



IV. STABILITY OF THE SATURATED STATE. 
PERIODIC DOMAIN STRUCTURE 

When the current overrides the value Jc the inclined 
uniform state (20) becomes unstable. The new stationary 
state, which occurs as a result of the instability, depends 
on the wire thickness h. If h exceeds some critical value 
he, which will be determined bellow, the nanowire very 
rapidly goes to the saturated regime where the magneti- 
zation is uniformly aligned along i-axis. This behavior 
is demonstrated in the Fig. 2b. It should be noted that 
contrary to infinite wires (or closed wires) the transition 
to saturation for wires of finite length is not quite a sharp 
jump due to the transitional formations of saturated do- 
mains at the wire ends. For thin wires with thickness 
h < he the saturated state is unstable with respect to a 
periodic domain structure formation, see region "11" in 
the Figs. 2, 3. Magnetization distribution of the domain 
structure is shown in details in the Fig. 4. Stability of 
the domain structure was proved using micromagnetic 
simulations. 

To determine the critical thickness he and the satu- 
ration current for h < he wc consider the stability of 
saturated state, when the quantization axis is directed 
along £-axis. Substituting = and <I> = 0^"^ into (15) 




FIG. 3. Diagram of stationary states of permalloy square 
nanowires with the transverse size h — hi under spin- 
polarized current influence, panel a). Sense of regions I, II, 
III and I-III is the same as in the Fig. 2. Panels b) and c) 
show the magnetization states for regions I and II respec- 
tively. Lines 1, 1', and 1" represent the critical current (25) 
for the cases A = 2, A = 1, and A — >■ oo respectively. Line 
2 shows the critical current (30). Crossing of the lines 1 and 
2 determines the critical thickness he whose dependence on 
parameter A is shown on the panel d), the point shows the pa- 
rameters of simulations: A = 2. Results of the micromagnetic 
simulations are shown by markers: disks and squares repre- 
sent the cases of the periodic boundary conditions and finite 
wire length L — 1 fim respectively; open and filled markers 
corresponds to transition from the inclined uniform state to 
the domain structure and to the saturated state respectively 
for the case of current increasing. Markers with dots shows 
the currents when the saturated state become instable for the 
case of the current decreasing. 

enable us to write the energy of the system in form 

S _ _ (26) 
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FIG. 4. The domain structure which appear in nanowires 
with transverse size h < he- The data obtained using micro- 
magnetic simulations. The 170 nm length part of a nanowire 
with thickness h — 5 nm and L — l^m is shown in the upper 
part of the figure. Magnetization distribution of the domain 
structure is shown by arrows and color scheme. Magnetiza- 
tion components along the wire are also plotted at the bottom: 
solid, dashed 
respectively. 



and equation of motion (9) in form 
. ; .r 7 



(27) 



where j = jT?/2. Since the Eq. (27) is linear we will look 
for solutions in form 



In this case the rate function C±(fc) is 
C±(fc) = -j±i^(A:), 

k{k) = ^j9g^k)~[e2k^~g{k)-2g{0)Y 



(28) 

(29a) 
(29b) 



As one can see from (29a) the current j plays the role 
of an effective damping as well as in case of the inclined 
solution (20). For j > >c{k) the rate function C± < 
and it means that the saturated state is stable. However, 
for j < maxfc i<(fc) the rate function becomes negative 
^+ < which means that the saturated state = is 
linearly unstable. Thus, the minimum current at which 
the saturated state is stable can be obtained as js = 
maxfc ii{k) or equivalently in physical units 



^ 8TTM^\e\h 

Js = max x(k) 

hr] k 



(30) 



Dependence of the critical current (30) on the transverse 
wire size h is shown by line 2 in the diagram of states. 



see Fig. 3a. The dependence Js{h) is not linear due to 
the nontriviality of dependence >c{h), see (29b) and (16). 

It is important to note that the appearance of the pe- 
riodic domain structure essentially depends on the re- 
sistance mismatch parameter A. Contrary to Jc the 
critical current Js does not depend on the parameter 
A. Using (16) and (29b), one can easily show that 

maxx(fc) — 3(7(0) = 1/4 for the case /i — > 0. Com- 
k 

paring now (30) and (25) we conclude that for A = 1 
(minimal possible value) the line Jc{h) is tangent to the 
curve Js{h) in point h = 0, see line 1' in the Fig. 3a. 
In this case Jc{h) > Js{h) and this means that the sat- 
uration, which appears when the current overrides the 
value Jc, remains stable and the domain structure is not 
formed. For the other case A > 1 the slope of the line 
Jc{h) decreases and the intersection of the dependencies 
Jcih) and Ja{h) determines the critical thickness he'- 



Jc{hc) = Js{hc) 



(31) 



such that for h < he the saturated state is unstable for 
currents J < Jg and the stable periodic domain structure 
appears in the interval Jc < J < Js- 

For h > he we have Jc{h) > Js{h) and a hystere- 
sis takes place. When current increases from zero and 
reaches the critical value Jc the inclined uniform state 
(20) becomes unstable and the system abruptly goes to 
the saturated state. By moving in opposite direction, 
i.e. when we start to decrease current from the saturated 
state it remains stable down to the current Js < Jc- Note 
in passing that the transition from the saturated state to 
the inclined state (20) may be accompanied by appear- 
ance of an irregular domain structure which arises due 
to degeneracy of the inclined state (20): if $ = $0 is 
a solution of (20) then $ = $0 + is another solution. 
The described hysteresis region is denoted "I-III" in the 
Figs. 2b), and 3a). 

Thus, under the adiabatic increase of the current the 
saturation value of the current depends on the wire thick- 
ness: 



Js 
Jr 



h < he, 
h > hr- 



(32) 



It is worth noting that for h > he the saturation current 
for wires of finite length is slightly lower than J^,, see 
Fig. 3a. This happens due to a particular role of the wire 
ends in the process of the saturation: for the currents J ^ 
Jc domains magnetized along £-axis occur at each of the 
ends. The length of this end domains rapidly increases 
with the current increasing and eventually the wire goes 
to the saturation state when the length of each of the 
domains reaches the half of the wire length. 

The critical thickness he determined by (31) depends 
on the parameter A, see Fig. 3d. With A increasing 
the quantity he asymptotically approaches its maximum 
value hf'''' « 2.303^, see hue 1" in the Fig. 3a. For 
the minimal possible value A = 1 the critical thickness 
/ic = 0. 
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V. CONCLUSIONS 

We study theoretically the influence of perpendicu- 
lar spin-polarized current on magnetization behavior of 
the long ferromagnetic nanowire of square cross-section. 
We found out that under the action of current the wire 
always reaches some stationary state. This is in con- 
trast to the case of planar films where either dynamical 
regimes or stationary states may realize for different cur- 
rent densities^^'^'^. In the no current case the wire is 
uniformly magnetized along its axis. For current increas- 
ing within the interval < J < Jc the wire magnetiza- 
tion remains spatially uniform being inclined in the plane 
perpendicular to the current direction. The inclination 
angle $ depends on the current density J. The critical 
current Jc depends on the thickness and the mismatch 
parameter A. When the current overrides the value Jc 
the nanowire goes either to the saturated state where 
the magnetization is uniformly aligned along the direc- 
tion of spin-polarization of the current or to a periodic 
multidomain state, depending on the wire thickness h. 
The domain structure exists for h < he in the interval of 
currents Jc < J < Js- In contrast to Jc the upper value 
of this interval Js does not depend on the resistance mis- 
match parameter A. 

For thick wires with h > he an opposite inequality 
Jc > Js holds and a hysteresis phenomenon takes place. 



where we omitted index n for the sake of simplicity. Solv- 
ing the system (Al) with respect to ip and -0* one obtains 



TO 

~D 



_9£_ 

dip* 



-ije TO 



dip* 



dip* 



(A2a) 
(A2b) 



Here D is the determinant of the system (Al): 



D = 



dm^ dm^ drri^ dm^ 



dip dip* dip* dip 



(A3) 



Substitution (5) into (A3) results in D = im^ and finally 
the equation (A2) takes the form of (6) with spin-torque 
term (7). 

Appendix B: Dipole-dipole contribution 



Substituting (5) into (14), we obtain that the harmonic 
part of the dipole-dipole energy in the wave- vector space 
in form of 



Appendix A: Equation of motion in terms of 
function 

In accordance to (4) one can consider that m„ = 
mn{ip, Ip*). Projecting the Eq. (1) to axes x and y re- 
spectively results in the system 

dm^ ■ dm^ ■ d£. ^ d£, . r z 

dip dip* dm^ dmy 



dmy ■ dmy ■ 9£ „ 9£ 

-Ip ~\- — — Ip — TO — — TO 



dip dip* 



dm^ 



dm^ 



£° - - 3v^sineao(0)4 + ^a,{k}M- 



-\ipk\^[a2{k) + 2a2{Q)]+c.c. 



where the coefficients ai{k) are 



ao{k) = 




EE 








ai{k) = 


a' 
SnJ'is 


EE 






&2{k) = 


SttJ^s 


EE 



cos Q {v? cos^ $ -I- sin^ $ — z^^) + i sin $ cos $ (n^ — y1 



ikn 
^ 1 



(n2 + + zl^) 



5/2 



vlii. - 22^^ - 3 sin^ e (n2 cos^ $ + yl^ sin^ $ - zl^) 



(n2 + yl^ + zl^) 



5/2 



(Bl) 



yl^) (cos 2$ + i cos e sin 2$) - sin^ 6 (n^ cos^ $ + yl^ sin^ $ - zl^) 



(B2) 



For the wires with cross-sections in form of square, disk 
or ring the indexes Uy and Vz are interchangeable, as well 
as indexes iiy and Hz- That enables us to to simplify the 



expressions (B2) 

Q!o(fc) = cos <I>(cos Q cos $ + i sin $)g(fc), 

ai(fc) = (cos 2$ - sin^ 9cos^ $ + i cos 9 sin 2$)5r(fc), 

a2{k) = (l-3sin^ecos2<I>).g(fc), 

(B3) 
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FIG. 5. Comparison of regular parts of sum (B4) for wires 
with square and round cross-sections. The solid and dashed 
lines correspond to the square and round wires respectively 
and represent expressions (B9) and (B12) respectively. 



where function g{k) is the following 



9 9 

^ — yf 



(B4) 



To obtain more suitable expression for function g{k) 
we proceed in (B4) from summation to integration. In 
this case g{k) can be presented as sum of regular and 
singular parts in the following way 



9{k) = greg{k) + gsing{k), (B5a) 
cxD h h 

9ren(k) = — — lim / dx cos(fca;) / dy / dzx (B5b) 

Tlh^ra^Qj J J 



X {h-y){h-z) 



7 7 



7r/2 a(x) 



gsing{k) = lim fdxcos{kx) [ dx [ drrx (B5c) 

Tlh^ro^Oj J J 



. A3- 



X {h — r sin x)(h ~ r cos x) 



2 2 2 

x — r COS X 

5/2 ■ 



Here we used the relation 

h h 



j dy j dy'F{\y - y'\) = 2 J dy(/i - y)F{y). (B6) 



Furthermore in (B5c) we proceed to the polar frame 
of reference (r, x) originated in one of the vertexes of 
the wire cross-section, and function cr(x) determines the 
cross-section shape. In the following we will ascertain 
that specific form of the function ct(x) does not affect 
the value of the integral (B5c). To do this one should 
change the variables r = rop, x = tq^ and find the limit 
in (B5c). The result is the following 

1 tt/2 oo 2 2 2 

gs^ng{k) ^ ^ J d^ J dx J <iPP ^ (^^2 ^ p2y/2 ' (^7) 




The direct integration of (B7) results in 

1 



9. 



Sing 



(fc) = - 



6' 



(B8) 



Let us now proceed to the calculation of the regular 
part (B5b). The direct integration over variables y and 
z with the next integration by parts over x enables us to 
simplify expression (B5b) 



greg{k) ^^^^ / sm{khx) 



(B9) 



V2^ 



1 



- 1 da; 



Using the representation (B9) one can easily obtain 
5reg(0) = 1/4 and consequently 

1 



5(0) = 



12' 



(BIO) 



Let us consider the case /ifc <C 1. In other words, the 
wire thickness is much smaller then the wavelength, it 
means that the form of the wire cross-section should not 
be of principle. To show this we calculate the regular 
part of function g{k) for a nanowire with disk-shaped 
cross-section with radius R = h/y/ir, so the cross-section 
areas of round and square wires are equal. In this case 
the regular part of the sum (B4) can be presented by 
integral 



2tv 2tt hi sJUhj s/tx 



9regik) 



47r/i2 

rr' cos{kx) [: 



ro 
X^- 



(Bll) 



(rcosx — r' cos x'Y 



[x^ + + 7-'2 — 2rr' cos(x — x')]^^^ 
Using the parametrization 



P 



1 



(p2-|-c2)5/2 (p2_,_c2)3/2 



OO 



and relation 

2ir 



Jo(v/^2 + r'2-2rr'cosx) = 27r Jo(r) Jo(r'), 



where Jq (x) is Bessel Function of the first kind, one can 
easily integrate the expression (Bll) and obtain 



9°regik) = ^^^l 



kh 

7^ 



kh 

7^ 



(B12) 



where Ii{x) and Ki{x) are modified Bessel functions of 
the first and second kinds respectively. Since the func- 
tions 9reg{k) and 9°f.g{k) are very close, see Fig. 5, one 
can use the approximation (16). 
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